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We present a new approach to gauge-invariant cosmological perturbations at second order, which 
is also covariant. We examine two cases in particular for a dust Friedman-Lemaitre-Robertson- 
Walker model of any curvature: we investigate gravity waves generated from clustering matter, that 
is, induced tensor modes from scalar modes; and we discuss the generation of density fluctuations 
induced by gravity waves - scalar modes from tensor perturbations. We derive a linear system of 
evolution equations for second-order gauge-invariant variables which characterise fully the induced 
modes of interest, with a source formed from variables quadratic in first-order quantities; these we 
transform into fully-fledged second-order gauge-invariant variables. Both the invariantly defined 
variables and the key evolution equations are considerably simpler than similar gauge-invariant re- 
sults derived by other methods. By finding analytical solutions, we demonstrate that non-linear 
effects can significantly amplify or dampen modes present in standard linearised cosmological per- 
turbation theory, thereby providing an important source of potential error in, and refinement of, the 
standard model. Moreover, these effects can dominate at late times, and on super-Hubble scales. 



I. INTRODUCTION 



The standard model of cosmology is built using perturbations of the homogeneous and isotropic Friedman- 
Lemaitre-Robertson- Walker (FLRW) models. These perturbations have, with a few exceptions, been limited to 
first-order thus far; at first-order, the generic large-scale features of the universe can be accounted for. However, 
there are good reasons for investigating higher-order perturbations as a refinement of the accuracy of this model. 
The unprecedented precision of present and upcoming cosmic microwave background (CMB) measurements, such 
as WMAP and Planck, are reaching levels where non-linear effects may play an important role, especially in any 
non-Gaussianity present. Indeed, it is rapidly becoming an important area of research to try to predict possible 
imprints of quantum gravity effects in the CMB; but this can only properly succeed if we also fully understand 
the range of classical non-linear effects which may produce very similar fluctuations. Non-linear effects may also 
play an important role in inflation as in the stochastic approach to inflation, which partly includes second-order 
back-reaction effects [1] , and also from a fully second-order approach [2] . 

Another reason to consider second-order perturbations is that there is no way, within linear theory, to decide 
when the perturbations have become too large for the theory to handle - i.e., when non-linearities should be taken 
into account. The reason is that, while it would seem sensible to compare the magnitudes of linearised objects to 
some invariant background variable, such as the energy density, say, this does not always work, a simple example 
being gravitational waves in flat space. The study of cosmology is not yet at an advanced enough stage for fully 
non-linear numerical predictions, as is the case for many astrophysical situations, so second-order perturbation 
theory is the tool which must be developed for this. 

For these reasons, second-order cosmological perturbation theory is beginning to be investigated, following work 
in the sixties by Tomita [3] who first investigated the second-order terms from scalar modes. For example, other 
authors have considered second-order perturbations of a fiat dust FLRW models, both with [4] and without [5] 
a cosmological constant. Inflation has also been considered at second order [2], providing the prediction of the 
bispectrum of perturbations from inflation. Meanwhile, the evolution of the curvature perturbation on super- 
Hubble scales after inflation has been shown to be constant at second order [6, 7]. More recently, it has been shown 
that second-order effects can lead to detectable non-Gaussianity in the CMB [8, 9], as well as polarization [10]. 
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In this paper we complement these results by demonstrating that second-order effects can dominate at late times 
and on super-Hubble scales over first-order modes. 

A crucial aspect of relativistic perturbation theory is the mapping, or gauge choice, between the (fictitious) 
background and the perturbed spacetime; many perturbation approaches are not invariant with respect to this 
mapping (see [11, 12] for a comprehensive discussion of this). In terms of coordinates, objects are typically not 
invariant under infinitesimal coordinate changes. This means that objects and modes specified in one gauge do 
not always correspond to objects and modes in another, making the interpretation of variables, and the physical 
situation, difficult. In fact, a gauge-dependent object is not observable [12], which follows from the principle of 
general covariance. It is therefore important to have a method of describing perturbations using variables which 
are invariant under gauge transformations. A gauge-invariant object at first-order is a quantity which vanishes 
(or is constant) in the background; this is the Stewart- Walker Lemma [11]. This gauge problem was first solved 
in cosmology by Bardeen [13] by providing a full set of gauge-invariant quantities with which to describe the 
perturbed spacetime. Bardeen's variables are formed from gauge-invariant (GI) linear combinations of gauge- 
dependent variables, given in a particular coordinate system; but this prescription can make their interpretation 
difficult. 

This situation has been clarified more recently by the covariant and gauge-invariant perturbation approach of 
Ellis, Bruni, and others [14, 15], using the covariant 1+3 approach [17]. Its strength in cosmological applications 
lies in the fact that it is well adapted to the system it is describing: all essential information can be captured 
in a set of (1+3) covariant variables (defined with respect to a preferred timelike observer congruence u a ), that 
have an immediate physical and geometrical significance. These variables satisfy a set of evolution and constraint 
equations, derived from Einstein's field equations, and the Bianchi and Ricci identities, which form a closed system 
of equations when an equation of state for the matter is chosen. The covariant and gauge-invariant linearisation 
procedure is easy and transparent: it consists of deciding which variables are 'first order' (or 'of order e') and those 
which are 'zeroth order' - i.e., those which do not vanish in the background, which is usually a FLRW model. 
Products of first-order quantities can then be ignored in the equations. The key point of the approach is that it 
deals with physically or geometrically relevant quantities, such as the fractional density gradient, T> a = aD a fi/ fi 
(where a is the scale factor, \i is the gauge-dependent energy density in the perturbed spacetime, and D a is the 
derivative operator in the observers' rest space), and the comoving spatial expansion gradient Z a = aD a 9 (where 
8 is the gauge-dependent expansion), being the most important for scalar perturbations, while the electric and 
magnetic parts of the Weyl tensor, E a b and H a b, respectively, represent the non-local parts of the gravitational 
field, and describe, amongst other things, the propagation of gravitational waves. All 3-vectors and projected, 
symmetric, trace-free tensors in this approach, are automatically GI because of the symmetry of the background. 

The gauge problem at second-order is quite a challenge: gauge transformations at second-order for second-order 
variables are huge in general - often as complicated as the evolution equations themselves (see e.g., [2, 5, 12]). 
This is a problem far greater than in linear theory in terms of extracting the physics of non-linear effects, simply 
because it is far harder to intuit the physics by considering a selection of different gauges; the equations are just 
too complicated, and the range of possible gauges vast. So even choosing a specific gauge at first and second 
order, one has no guarantee that one has fully understood the problem. The key problem is that for an object to 
be gauge-invariant at second-order it must vanish in the background and at first-order (or be constant) [12]. This 
problem has hints of a solution in the recent work of Nakamura [18], using the methods given in [12], in which a 
general formalism for deriving GI quantities at second-order, assuming the existence of a method at first-order, 
is given. Thus it is possible that a fully gauge-invariant second-order theory may be presented in the near future 
using Bardeen's method for constructing first-order GI quantities together with Nakamura's approach. There has 
been recent progress along these lines; [2] have given a gauge-invariant definition of the comoving curvature at 
second-order for scalar perturbations. Nevertheless, it is likely that when such a theory is found it will be very 
difficult to understand exactly what everything means, especially as the construction is not unique, although it 
will likely be crucial for numerical prediction of observable quantities such as the CMB power spectrum. 

In this paper we discuss an alternative route to second-order gauge-invariant perturbations in cosmology 
following the covariant route initiated in [14]. The covariant 1+3 approach is a good method for perturbation 
theory simply because it is manifestly covariant; many of the usual gauge-problems arise as a result of coordinate 
transformations at each order. Once we have a complete set of GI variables, then we know that they are 
observable [12], and correspond to physical or geometrical quantities. 1 



1 It should be noted that there is a freedom inherent in the 1+3 approach, even after a complete set of GI variables have been specified, 
and this freedom lies in ones choice of observers at the relevant perturbative order: but this is not a gauge choice in perturbation 
theory, as it has nothing to do with ones mapping between the background and perturbed spacetime; it is simply a choice of 
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Using a zeroth-order background of a dust FLRW model with zero cosmological constant, we discuss 'mode- 
crossing' perturbations; that is, perturbations at second order generated from an 'orthogonal mode' at first order - 
second-order tensor modes generated by first-order scalar modes (Section III A), and viceversa (Section IIIB). 
2 We define a complete set of second-order gauge- invariant quantities which fully describe the induced modes, 
and convert all products of first-order quantities appearing in the second-order equations into new second-order 
variables, thus making the equations explicitly GI at second order, in a manner very similar to [16]. This 
simplifies the presentation of the equations by converting them into a linear system of DE's, a procedure which 
also simplifies their solution. In the case of first-order tensor modes, this system is infinite dimensional reflecting 
the coupling of modes of each wavelength feeding into the induced scalar modes. The standard zeroth-order 
harmonic functions can be used on second-order variables to remove the tensorial nature of the equations, and 
remove spatial gradients so converting the system into odes which may then be easily integrated. We explicitly 
integrate both cases under investigation to show that scalar-mode coupling effects may play an important role at 
both early and late times, and on small and large scales. 



II. FIRST-ORDER PERTURBATION: THE COVARIANT FORMULATION 



The 1+3 approach relies upon the introduction of a family of observers travelling on a four- velocity u a , with 
which all geometrical and physical objects and operators - essentially the Riemann curvature tensor and the 
covariant derivative - are decomposed into invariant parts; scalars along u a and scalars, 3- vectors, and projected, 
symmetric and trace- free (PSTF) tensors orthogonal to u a , as well as an evolution derivative along u a and a 
spatial derivative orthogonal to it. The Einstein field equations are supplemented by the Ricci identities for u a 
and the Bianchi identities, forming an complete set of first-order differential equations. We refer to [17] for details 
and references; we also follow their notation and sign conventions 3 . 

Perturbations at first-order around a FLRW background are relatively straightforward: define the first-order 
gauge-invariant (FOGI) variables corresponding to the spatial fluctuations in the energy density and expansion: 

X a = D aM , (1) 
z a = B a 6. (2) 

These are FOGI simply because they vanish in the exact FLRW background [11, 12]. We shall not use the 
more physically motivated variables V a = aD a ^/n and Z a — aD a 9, as these tend to make the equations more 
complicated to derive at second-order. Ignoring rotational perturbations the set of FOGI equations consists of 
the evolution equations 

(3) 
(4) 
(5) 
(6) 

(7) 
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- curl Hab = 
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+ curl Eab = 


— OHab', 





observers in the perturbed spacetimc, a choice one must make in any non-vacuum spacetime. For example, the gauge-invariant 
density gradient at first-order, T> a , can be set to zero by a first-order Lorentz boost, at the expense of introducing acceleration, 
say. Under a zeroth-order boost, however, V a would no longer be gauge-invariant as it would not vanish in the background, and 
would correspond to unusual observers in the background. In metric-based approaches (where the components of the metric are 
explicitly solved for), one must specify a perturbed velocity relative to the background, as well as having this frame freedom. This 
is because the Stewart- Walker lemma and its generalisations refer to the principle of general covariance, and not the principle of 
relativity: different observers measure different values for physical quantities such as the energy density, the magnetic field or the 
electric part of the Weyl tensor, for example. The only objects which are observer (or tetrad) independent arc the metric, the 
Riemann tensor (and its trace and trace-free parts) and the Maxwell tensor. 

2 We ignore rotational modes throughout. Only if they are included at first-order will rotational modes occur which do not occur at 
first-order, in the absence of acceleration. 

3 We use the standard notation whereby a dot represents differentiation along the observers' four-velocity u a , ip a ---b = u, c S7 c ip a ...b, 
and D a is a derivative in the rest space of the observers; D c -!/> a -.-6 = h c d h^ ■ ■ ■ h. V^i/V--/; where h a f, = g a ^ + u a M is the usual 
projection tensor orthogonal to u a . We use angled brackets on indices to donate the projected, symmetric and trace-free part of 
a tensor. We define the divergence of a vector t(> a as div ip = T) a tp a , and of a PSTF tensor %p ab as div i/> a = T> b ip ab ; and we define 
the curl of a vector as curl ipa = e a i yc D h 'ip c , and of a PSTF tensor, curl i/>af> = e cd ^ a D c ^ 6 ^ d ; e aoc is the observers' rest-space volume 
clement. 
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and the constraints 



= 


\X a - d\vE a , 


(8) 


= 


\Z a - divcr a , 


(9) 


= 


diviJ a , 


(10) 


= 


H ab - curl a ab , 


(11) 


= 


curl X a , 


(12) 


= 


curl Z a . 


(13) 



The whole system of equations is governed entirely by the shear, which obeys the covariant wave-like equation: 

a ab + curl curl a ab + ^6& ab + (|6» 2 - §/x) a ab = 0. (14) 

Each other object is determined by solutions of this equation - no further integration of the equations is required 
after the solution for the shear is found. Indeed it is useful to use {<7 ab ,& ab } as the 'basis vector' for the full 
solution, as we shall do later. While each variable in the problem may be shown to satisfy a wave-like equation 
of some sort, these do not necessarily close; in particular, the wave equation for E ab in the case of gravity wave 
propagation has a source term from the shear [19, 21]. 

Because the background is homogeneous and isotropic, each FOGI vector may be uniquely split into a curl-free 
and divergence-free part, usually referred to as scalar and vector parts respectively, which we write as 

V a = V Sa + V Va where curl Vs a = and div V v = 0. (15) 

Similarly, any tensor may be invariantly split into scalar, vector and tensor parts: 

Tab = T Sab + T Vab + T Tab where curlTs Qb = 0, divdivT v = 0, and divT Ta = 0. (16) 

In the equations we can separately equate scalar, vector and tensor parts. As we are ignoring rotation, all vector 
parts are zero. 

Each FOGI vector consists of only the scalar part and is curl-free; so 

X a = Xs a , (17) 

similarly for Z a . Each tensor consists of a scalar and a tensor part, which describe clustering matter and gravity 
waves respectively. Thus the equations can still be manipulated further: splitting the shear into its curl free and 
divergence free parts gives us an oscillatory equation for the scalar part, which tells us the increase in shear found 
as matter clusters: 

<*« 6 + §0A o6 + (le 2 - |m) K ab = o, (is) 

and a wave equation for the tensor part, which is the equation governing gravity waves in an expanding dust 
universe: 

*a6 - D'^afc + f 0*„6 + + | M ) *r o6 - 0, (19) 

which may be shown using the identity for a FOGI tensor: 

curl curl T ab = -D 2 T ab + §D (a div T b) + ( M - \d 2 ) T ab . (20) 

The other two tensors can be similarly split. Because the magnetic Weyl curvature is the curl of the shear, the 
scalar part of it vanishes at this order. 

We shall use Eqs. (18) and (19) to form the source for second-order perturbations. 

III. SECOND-ORDER PERTURBATION 

Lets consider non-linear perturbations up to second order in the shear, ignoring rotation and fluid modes at 
both first and second order, but keeping both scalar and tensor modes. The shear 'wave' equation becomes 



a ab + curl curl <j ab + \6a ab + (p 2 - a ab = o c{a o h) c + 9<r c{a a b) c 7 



(21) 
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where the 'curl-curl' operator now has the quadratic contributions when converting to the Laplacian: 

curl curl T ab = -V 2 T ab + §D (a div T b) + (/x - \0 2 ) T ab + 3E c{a T b) c - 8a c{a T b) c . (22) 

In deriving Eq. (21) we have substituted for the Weyl curvature in terms of the shear, given by Eqs. (5) and (11). 
We have neglected terms 0(a 3 ) (and derivatives of), which is a consistent covariant linearisation procedure to 
second order. There are some important things to note about Eq. (21): 

1. It is not actually second order in the usual perturbative sense. The shear appearing in the quadratic shear 
terms on the rhs should not be the same as that on the left; on the left we have a mixture of first and second- 
order quantities, but on the right, the terms are made up of products of first-order variables. Integrating 
Eq. (21) without explicitly taking into account this distinction will give a different solution than explicitly 
making it so [by using solutions of Eq. (14) for the terms on the right]. 

2. It is not gauge-invariant. The mixture of first and second-order quantities on the left is really the core of the 
gauge problem: the 'first-order bits' don't cancel out upon application of Eq. (14) because the derivative 
operators in each of the two equations are not the same. In the second-order equation there are derivative 
operators of order one hanging around; but solution of the linear problem does not tell us what these are. 
Hence, the equation can't be integrated. 

We shall now convert Eq. (21) into a gauge-invariant equation at second order in the two cases we have 
mentioned. 



A . GRAVITY WAVES FROM DENSITY FL UCTUATIONS 



Consider the case where we excite only the scalar modes at linear order, which describe density perturbations. 
These modes are covariantly characterised at linear order by the GI condition 

curl <j ab = 0. (23) 

Hence, we may define the variable 

S ab = curl <j ab (24) 

which is second-order and gauge-invariant up to and including second-order (SOGI), as it vanishes at all lower 
orders [12]; this variable forms the core of the analysis of gravity waves generated by density fluctuations. In fact, 
it is just the magnetic part of the Weyl tensor if rotational modes are ignored (as may be seen from the curl a ab 
constraint); nevertheless, it is a distinction worth keeping, as this would not usually be the case. Taking the curl 
of Eq. (21), or taking the time derivative of the magnetic Weyl evolution equation results in the equation for T, ab : 

£ afc + curl curl S Qb + |(9£ afc + ((9 2 - fi) S a6 = S ab (25) 

where the source term is given by 

S ab = e cd{a {-\a h) d d\vo- c - |cr b) d div cr c - pa b) d d\v<7 c 

• ec I ^ ^ d o^- ec T\ ' d I 1 „dey\c • i 1 •dej^C' 
— a D| e |Cfy — lO D| e |0f,) + 2 a U a b)e + ^ U ab ) e 

+ ia 6) e DV + f <7 6> e DV + 0D C {a b)e a de )}. (26) 

This source term is considerably more untidy than that of the gauge-dependent equation for the shear, Eq. (21) 
because of the extra terms which arise from successive applications of the commutation relation 

(curlT afc )' -curlT afc = -i0curlT ab -e cd<a( T ec D| e |T fc) d , (27) 

which holds for this case. 

The second-order equation (25) is not yet properly SOGI because the rhs contains terms which are FOGI, 
although the rhs as a whole is second-order, and S ab is a SOGI tensor. To make the source explicitly composed 
of SOGI quantities, we define the SOGI variables: 

ipiab = e cd(a (divcr c )cr fc) d , (28) 

foab = e C(i <a(divcr c )(T fc) d , (29) 

tpsab = e cd{a (d\va c )a b) d , (30) 

"04ah = e cd(a (divo- c )o- b) d ;. (31) 
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It is straightforward to show that, to second order, these variables satisfy the closed system of evolution equations 

^ab = *V> 6, (32) 
where ip ab = (V>ia6, V^ab, tp3ab, i^Aab), and the scalar-scalar coupling matrix is given by 



/ -\6 1 1 \ 

-!<9 2 + | M -26 1 

-§0 2 + §/x -26 1 



(33) 



which is derived using Eq. (18). Similarly, we may define the variables 

4>iab = e cd (a{d\so- c b)e )a de , (34) 

<t>2ab = t cd{a {A\S(T C b)e )cT de , (35) 

foab = e cd{a (d\sa c b)e )a de , (36) 
4>4ab = e cd{a (d\s& c b)e )a de , (37) 

where we have defined the distortion of a ab as 

dis <T cah = D( c er a6 ), (38) 

which is the remaining invariantly defined part of the spatial derivative of a PSTF tensor which is not part of 
the divergence or curl [22]. It is clear that these variables, <f> ab , obey the same evolution equation as ip ab , given 
by Eq. (32). By performing a complete 1+3 split of the derivatives of the shear in the source of Eq. (25) to 
write them in terms of the divergence and distortion only, we may write our source term in terms of these SOGI 
variables: 

Sab = -^tylab ~ %i>2ab + ^i>3ab + 0<t>lab + %<hab + %<l>3ab- (39) 

We now have a complete closed linear system of differential equations which may be expanded in the usual 
background harmonics, and integrated straightforwardly (especially numerically) once initial data is prescribed. 
By taking various spatial gradients of these eight variables it can be shown that there are some (differential) 
constraints which must be satisfied amongst these variables, which serve to constrain initial data, but we shall 
not pursue this further here (see [23] for details of where these come from). Note that we have managed to deal 
with all the complicated scalar-scalar coupling without recourse to harmonics. 

The mode-mode coupling variables, ip ab , <$> ab do several things in aiding the solution to Eq (21) at second-order: 
they transform Eq. (25) into a genuinely second-order equation in the perturbativc sense because in calculating 
their evolution equations we neglected terms 0(a 3 ), and made the final distinction between the variables on the 
lhs and those on the right; they turn Eq. (25) into an explicitly SOGI equation consisting only of SOGI variables; 
and finally they turn the system of equations which must be solved into a system of linear differential equations, 
which are much easier to integrate numerically, and analyse by other techniques, than systems which are not 
linear. 



1. THE INDUCED TENSOR MODES 



The mode coupling which provides the source for the second order perturbation, although being products of 
scalar mode variables, induces tensor modes at second order (and, in more general situations, rotational modes), 
seen by the fact that curl S ab ^ 0. This mode-mixing is perhaps one of the most interesting parts of non-linear 
perturbation theory; as we show here, for example, gravity waves may be generated - or reduced - purely from 
the effects of linear clustering of matter. The second-order scalar modes which T, ab tells us about, however, are 
not genuine second order effects: for scalar modes, divS a ^ 0, and is in fact just a combination of quadratic 
first-order variables, as may be seen through the div_ff a constraint. It is interesting to note that 

divS a = e abc a bd a c d = ^ <T(k)V(k') e abcQ\t)Q (k') C d , (40) 

k,k' " v ' 

=o when k=k' 
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using standard harmonics [17, 19], so coupling between modes of differing wavelength do not give rise to tensor 
modes at second order, for which div£ a = = div H a . This is a reflection of the integrability conditions 
discussed in [23], which suggests that div H a = may not be as restrictive in perturbative models as it is in the 
fully non-linear case. 

The second order tensor modes, on the other hand, are covariantly characterised by div S a = 0, for which we 
have the genuine wave equation: 

% ab D 2 ^ ab + l6% ab + I9 2 ^ ab = S Tab , (41) 

upon using Eq. (22) in Eq. (25). Harmonic functions, defined as usual in the background (any changes would 
be 0(3)), may now be introduced, and the equations integrated in time; we discuss this in the case of a flat 
background below. 

Now, gravity waves are waves in the Weyl curvature which induce waves in the shear (non-local waves which 
travel at the speed of light and thus can't be removed by a change of observers) , so in order to fully understand 
gravitational radiation induced by density fluctuations, we must relate £ a & to the Weyl curvature. This simply 
involves creating SOGI variables for the Weyl curvature. The magnetic part is easy, as we have H ab = T, ab in the 
absence of rotation, while for the electric part of the Weyl tensor we may define 

S ab = curl E ab = -± ab - 0Z ab + f Viafc, (42) 

which completes the solution for the induced tensor modes. 

In the case of a flat background, the equations are particularly easy to integrate (see [19] for details on solving 
the homogeneous part of Eq. (41) in a non-flat background): we have the scale factor a = dot 2 / 3 (we set a = 1), 
implying 9 = 2/t and /z = |# 2 ; integrating Eq. (32) for xj) ab and (p ab , we find that the source has three components: 

Sab = a ab t- J / 3 + f3 ab t- 4 + 7 aht- 17/3 , (43) 

where the last term is the largest at early times (which has the time dependence of a a&). The tensors a a b, Pab 
and 7 a b are constant in time and are a combination of the initial conditions of ip ab and <fi ab ; the details of this 
need not concern us here. Using standard tensor harmonics [15, 19] removes the tensorial nature and spatial 
dependence of the wave equation (41), and replaces the covariant Laplacian with a harmonic index, D 2 — k 2 /a 2 , 
the full solution being a sum over harmonic modes. The analytic solution to Eq. (41) may be easily obtained: 

= £< fc) + + a^k-H- 1 + p^k-H-*/ 3 + 7 (fc) x mess (44) 
where the two solutions to the homogeneous part are 

E (fc) = s W t -a/3 ^ 3x 2 _ -q cos3a , _ 3xsin3:z .] (45) 

E< fc) = S< fe) i- 8/3 [(3a; 2 - l) sin 3x + 3x cos 3a;] (46) 

where x = kt 1 ! 3 and s[ k \ are constants. The 'mess' from the j a b term dominates the solution at early times 
and on super-Hubble scales, when tel and x <C 1 respectively: 

E« ~ | 7 (fc) t - n /3 [1 - fx] + 0(t- 8 / 3 ), (47) 

with the next term arising from the homogeneous solution. Meanwhile, at late times we have a new power law 
behaviour dominating from the scalar mode interaction 

£W ~ fc-V^- 1 + 3k 2 t- 2 {S { 2 k) sin3x - + i^A; 3 7r 7 (fe) ) cos3a;} + 0(t- 7 / 3 ), (48) 

while we see that on small scales (k ^> 1) at late times the oscillatory behaviour is also dominated by the scalar 
mode interaction. Thus, in many regions of interest the terms arising from the scalar mode coupling dominate 
the solution. 

In Fig. 1 we show some solution curves for T,^ for a flat model both with and without the non-linear source. 
This shows that the effect of including the scalar-scalar interaction in the generation of gravity waves may increase 
or decrease the amplitude and power of the induced tensor perturbations for a given wavelength, possibly quite 
significantly. This may have implications for the spectrum of tensor modes in the cosmic microwave background 
(see e.g., [19, 20] for details of calculating the CMB power spectrum from the solutions given here). 
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FIG. 1: The induced gravity waves from non-linear density fluctuations for a flat background. Here we have shown the 
harmonic component of H a b, T,^ k \ corresponding to k/ao = 25 in three situations, all of which have initial conditions such 
that E^ = 1 and E' fc ' = at t = 0.1 (in arbitrary units). The thick grey curve represents the oscillatory decay of gravity 
waves with no source present, as predicted from linear theory. To contrast with that we have shown two different cases: 



curve Si has parameters from the source term a' fe ) = 0, /3' fc ' = 0.8, 
of the gravity waves over the linear case, while curve S2 has a' fc ' 
damping. 



„(*0 



= —.03, which shows considerable amplification 
(5 {k) = 0, 7 (ft) = .02, which shows significant 



B . DENSITY FL UCTUATIONS FROM GRAVITY WAVES 



We shall now look at the reverse situation from the previous section, and consider the case where only tensor 
perturbations are excited at first order. To second order we shall still neglect rotation and fluid modes for 
simplicity of exposition, and study the induced density fluctuations, characterised by the scalar modes. 

Pure tensor modes are covariantly characterised by the FOGI condition 

divcr a = (49) 

from which it follows that the divergence of the two Weyl curvature tensors must vanish also. Thus, at second 
order, we use the SOGI variable 

A Q = div(7 a (50) 

as our key variable. From the div i? a -constraint and the non-linear commutation relation 

div curl T a -\ curl div T a = e abc E bd T c d - \9t abc (r hd T c dl (51) 

we find that curl A a = because we're ignoring rotation - so A a contains only scalar modes. Our wave-like 
equation for a ab , Eq. (25), becomes, upon taking its divergence, 



A a + l9A a + (6 2 - A a = T a 



(52) 
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where the source arising from tensor-tensor coupling is, neglecting terms 0(a 3 ) and 0(a x A), 



% = -§0D O (a bc a bc ) - |D a (a fcc( x bc ) , 
which is remarkably simple, despite several applications of the commutation relation Eq. (51) and 



(53) 



(divT„)" - divT a = -l6d\vT a -a bc D b T ac + e abc H bd T c d . 



(54) 



We have again substituted for the Weyl curvature in terms of the shear. 

As in the previous case, we may make the source for Eq. (52) explicitly SOGI by introducing a set of SOGI 
variables; variables which are constructed from quadratic FOGI variables. However, the presence of the Laplacian 
in the wave equation for the tensor modes at first-order means that we need an infinite set of these variables; a 
mode of one particular wavelength at first-order couples with all higher wavelength modes at second order. One 
can approach this in two ways: explicitly split the first-order solution into a sum over tensor harmonics; when 

inserted into the source term T a we have an infinite sum (over k and k') over terms of the form a^a^ k ^Q^Q^ k \ 
each of which may be expanded as an infinite sum over tensor harmonics. Or, one can covariantly define an infinite 
set of SOGI quadratic variables, without expanding the first-order solution in harmonics. Both methods are 
more or less the same when the tensor harmonics form a complete set of basis functions over which any (smooth) 
solution of the first-order wave equation (19) may be expanded. However, there are certain situations in which 
the first-order solution may not be completely given as a sum over harmonic modes - black hole perturbation 
theory being perhaps the most widely known case [25] 4 . Therefore we will explore the second route, although 
the first may be more suitable for many applications. 
We define the variables 



1 /,( m '") 

Via 



1 /,( m ^™) - - 

r3a 



_ 2n+2m+li 



a ' L D 

„2n+2m+l r > 
a 2n+2m+l D 



Via 



D( 2ro V 6c D( 2 "V fec 
D( 2m V 6c D( 2 ")<7 6c 



= a 2 "+ 2m+1 D a D^a bc D^a bc 



(55) 



where the indices m,n are positive integers (including zero), and D^ 2 ™) represents the Laplacian operating n- 

times. Note that ipi and ^4 are symmetric on the m and n indices, while ip2a' n ^ = V4a'"^ ~ these forming a set 
of constraints. The variables evolve, to second order, as 



(56) 



where the tensor-tensor coupling matrix is given by 



/ 

l/)2 1 





1 

-f* 



9^ - 6^ 



1/J2 



\ 
1 

1 



(57) 



while the cascading matrices are 



A = 



/ : 2 • \ 

a 2 ■ 



V 



• a- 2 • J 



X = 



■ J 



(58) 



which couple modes of differing wavelength and feed them into the source for the scalar modes. We also require 
the set of variables 



,(m,n) 
'\a 



a 



2n+2m+ljy b 



D( 2m V c(a D( 2 "V h) c 



4 Expanding the key wave equation for first-order gravity waves around a black hole in temporal harmonics does not sum to the 
complete solution, because of purely outgoing boundary conditions in that problem. 
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^(m,n) = a 2 n+ 2 m +l D 6 
(m.n) _ a 2n+2m+lY)b 



(m,n) _ /i 2n+2m+l-|- ) b 



u 4a 



= a 



D( 2 ™)<7 c<a D( 2 ")<7 fe) c 
DMa c(a DWa, 



6) 



(59) 



where is is clear that these variables (j)^ n,n ^ obey exactly the same evolution and constraint equations as ip < £ n,n \ 
Our source now becomes explicitly SOGI, so that the full evolution equation for induced scalar perturbations 
from gravity waves becomes very simple: 



A a + l6A a + (6 2 - M ) A a = -§0a" Vi„ ! 



(0,0) 



(0,0) 



(60) 



We now have a SOGI linear system of pure evolution equations, which may be integrated in a straightforward 
manner (although the infinite-dimensionality of the system may pose a few problems; this would typically be 
dealt with with a short wavelength cutoff) . There will also be a system of constraint equations between all the 
</>i m ' n) and Vi 7 "'^ which are derivable from their definition; these would serve to constrain initial data when 
integrating the full system of equations (or, if the system of equations were analysed in tensor harmonics, would 
provide a system of algebraic equations which could then be solved, giving a smaller system of evolution equations 
to integrate). 

The complete solution for scalar modes is given when Eq. (60) is integrated: the spatial gradient of the 
expansion is given by the div-shear constraint 



Z a = 3 A a 



while the gradient of the energy density is a little more complicated: 

X a = -3A a - 36A a + 2a-Vi°„' 0) + a" Vi°' 0) 



(61) 



(62) 



follows from the div E a constraint. 

We may once again integrate our key equation (60) in a flat background. As an illustration, consider only 
super-Hubble tensor modes, and set m = n = 0. The source term may be found by integrating Eq. (56), and has 
the same time dependence as the source in the previous case, Eq. (43). The solution for A a is 



A a = T la r 1 + T 2a t-«/ 3 + ^a a r 1 / 3 - \p a t 



7a t 



-11/3 



(63) 



INDUCED MODES 



which shows that the tensor mode coupling dominates at both early and late times, and, in fact, contains a 
growing mode (A a /9 ~ aa). 

At late times density perturbations are the key for structure formation in the observable universe, and a 
relativistic analysis must be used for scales approaching the Hubble length. Density perturbations may be directly 
related to the growth rate of clustering matter [26] . As tensor modes may alter the behaviour of clustering matter, 
evidence may be found in variables such as the bias parameter. 



IV. CONCLUSIONS 



Understanding non-linear effects in cosmology will likely play a pivotal role in cosmology in coming years, as 
observations pin down the overall properties of the universe ever more tightly. Most of the problems preventing 
a systematic study of non-linear relativistic effects are mathematical: fully non-linear relativistic simulations are 
beyond the horizon in cosmology at present, leaving perturbation theory as our main probe. But non-linear 
perturbation theory is not particulary simple as gauge problems and large equations make it quite untidy. In this 
paper we have considered for the first time the 1+3 covariant approach to gauge-invariant non-linear perturbation 
theory 

The main problem the 1+3 covariant approach has at second order is intricately related to the gauge-problem: 
while second-order equations, such as Eq (21), can be written down - by crossing off all the third order terms - 
they can't in general be integrated. This is because the derivative operators, through the commutation relations, 
contain lots of first order stuff, which the covariant approach does not solve for at first-order (it only solves for 
physical variables, not operators, in contrast with the metric approach). In order to integrate the second-order 
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equations therefore, derivative operators (dot and D) must operate only on variables which vanish at first and 
zeroth order - SOGI variables, in other words. We created these SOGI variables in two ways: by appropriate 
derivatives of first-order quantities which vanish at first-order - exactly the same trick used to create FOGI 
variables from zeroth-order scalars; and from explicit products of FOGI quantities. 

The GI approach we develop here has the advantage of producing a relatively simple linear system of odes, 
which is easy to integrate. By integrating the SOGI equations, we have found that non-linear effects may well 
play an important role at both early and late times. Indeed in both cases we examined we found that second-order 
modes will dominate the perturbation spectrum at late times, which is a consequence of the fact that second- 
order modes arise partly as an integrated effect, and can't always be neglected by assumption. In addition, tensor 
modes generated by non-linear effects may play an key role on super-Hubble scales. This analysis suggests that 
linear perturbation theory may not be a sufficiently accurate tool to understand some subtle dynamical aspects 
of the universe, but further investigation is required to find out exactly what these are. 
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